Strongest gravitational waves from neutrino oscillations at supernova core bounce 
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Resonant active-to-active (v a — > v a ), as well as active-to-sterile (v a — > v s ) neutrino (v) oscilla- 
tions can take place during the core bounce of a supernova collapse. Besides, over this phase, weak 
magnetism increases antineutrino (v) mean free paths, and thus its luminosity. Because the oscilla- 
tion feeds mass-energy into the target v species, the large mass-squared difference between species 
(v a -4 i/,) implies a huge amount of energy to be given off as gravitational waves (^GWs ~ er S 
s _1 ), due to anisotropic but coherent v flow over the oscillation length. This asymmetric i/-fiux is 
driven by both the spin-magnetic and the universal spin-rotation coupling. The novel contribution 
of this paper stems from 1) the new computation of the anisotropy parameter a ~ 0.1 — 0.01, and 2) 
the use of the tight constraints from neutrino experiments as SNO and KamLAND, and the cosmic 
probe WMAP, to compute the gravitational-wave emission during neutrino oscillations in supernovae 
core collapse and bounce. We show that the mass of the sterile neutrino v s that can be resonantly 
produced during the flavor conversions makes it a good candidate for dark matter as suggested by 
Fuller et al. (2003). The new spacetime strain thus estimated is still several orders of magnitude 
larger than those from v difussion (convection and cooling) or quadrupole moments of neutron star 
matter. This new feature turns these bursts the more promissing supernova gravitational-wave 
signal that may be detected by observatories as L1GO, VIRGO, etc., for distances far out to the 
VIRGO cluster of galaxies. 



I. INTRODUCTION 

Supernovae neutrinos and gravity waves. — That outflowing neutrinos (vs) from a supernova (SN) generate gravita- 
tional waves (GWs) was firstly pointed out by Epstein (1978). However, over the first ~ 10 milliseconds (ms) (Mayle, 
Wilson & Schramm 1987; Walker & Schramm 1987) after the SN core bounce the central density gets so high that 
no radiation nor even vs can escape, they are thus frozen- in and strongly coupled to the neutron matter (N ) as 
described by the Lagrangean (see Kusenko & Postma 2002 for this dynamics) 



= ^ [JV° 7m (1 - 75)^V°] {^7 M (1 " 75)^} 
with the v field (tp) satisfying the time-dependent Dirac equation 
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At this phase the whole proto-neutron star (PNS) dynamics is dominated by gravity alone, and can be appropriately 
described by the general relativistic Oppenheimer-Volkoff equation for both the N° + v fluid (see Mosquera Cuesta 
2002). As discussed by Mayle, Wilson & Schramm (1987); and Walker & Schramm (1987), it is over this early transient 
that most v flavor conversions are expected to resonantly take place and consequently the super strong GWs burst 
from the oscillation process to be released. GWs from this decoupling has been suggested to likely be the ultimate 
process responsible for the neat kick given to a nascent pulsar during the SN collapse (Mosquera Cuesta 2000; 2002). 

The contention of this paper is a) to pave, in the framework of general relativity (GR), the pathway to this 
fundamental astrophysical process of generation of GWs from v oscillations in a PNS. b) to demonstrate, by taking 
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into account experimental and observational constraints, that v oscillations during SN core bounce do produce GWs 
of the sort predicted by Einstein's GR theory, and more crucial yet, c) to stress that these bursts are the more likely 
SN GWs-signals to be detected by interferometric observatories as LIGO, VIRGO, GEO-600, etc. We speculate that 
such a signal perhaps might have been detected during the SN1987a event, despite the low sensitivity of the detectors 
at the time. Some claims in this direction were presented by Aglietta, Amaldi, Pizzella, et al. [2], and related papers. 

II. THE MECHANISM FOR GENERATING GWS DURING NEUTRINO OSCILLATIONS 

To start with, let us recall how the production of GWs during v oscillations proceeds by considering the case of 
oscillations between active and sterile neutrinos in the supernova core. The essential point here is that oscillations 
into sterile neutrinos change dramatically the energy and momentum (linear and angular) configuration of the system: 
neutrinos plus neutron matter inside the PNS (check Eq.(2)). In particular, flavor conversions into sterile neutrinos 
drive a large mass and energy loss from the PNS because once they are produced they freely escape from the star. 
The reason: they do not interact with any ordinary matter around, i.e., they do couple to active v species but neither 
to Z° nor to vector bosons. This means that oscillations into steriles, in dense matter, take place over longer 
oscillation lengths, compared to v a — > v ai and the steriles encounter infinite mean free paths thereafter. Physically, 
the potential, V s (x), for sterile neutrinos in dense matter is zero. In addition, their probability of reconversion, still 
inside the star, into active species is quite small (see discussion below). This outflow translates into a noticeable 
modification of the PNS mass and energy quadrupole distribution, which as discussed below is dominated from the 
very beginning by rotational and magnetic field effects. 

Since most steriles neutrinos escape along the directions defined by the dipole field and angular momentum vectors 
(see Fig. 2), then the v outflow is at least quadrupolar in nature. This produces a super strong gravitational-wave 
burst once the flavor conversions take place, the energy of which stems from the energy and momentum of the total 
number of neutrinos participating in the oscillation process 1 . Further, the gravitational- wave signal generated this 
way must exhibit a waveform with a Christodolou's memory (Mosquera Cuesta 2002). 

The remaining configuration of the star must also reflect this loss. Hence, its own matter and energy distribution 
becomes also quadrupolar. Because this quadrupole configuration (the matter and energy still trapped inside the 
just-born neutron star) keeps changing over the time scale for which most of the oscillations take place, then GWs 
must be emitted from the star over that transient. At the end, the probability of conversion and the v flux anisotropy 
parameter (a, see below) determine both how much energy partakes in the process and the degree of asymmetry 
during the emission. Both characteristics are determine next. 

EDITOR PLACE FIGURE 1 AROUND HERE !!!!! 

The case for oscillations among active species is a bit different, the key feature being that mass and energy is 
relocated from one region to another inside the PNS, especially because of the weak magnetism of antineutrinos that 
allows them to have larger mean free paths (and thus oscillation lengths) (Horowitz 2002). In addition, oscillations 
of electron neutrinos into muon or tauon neutrinos leave these last species outside their own neutrinospheres, and 
hence they are in principle free to stream away. These neutral-current interacting v species must be the very first 
constituents of the v burst from any supernova since most v e s are essentially trapped. This must also generate 
GWs during that sort of flavor conversions, although their specific strength (strain) must be a bit lower compared to 
conversions into sterile neutrinos where almost all the v species may participate, and the large Am 2 in the process. 

Since the sterile neutrinos escape the core over a time scale of a few ms, the number of neutrinos escaping and their 
angular distribution is sensitive to the instantaneous distribution of neutrino production sites. Since thermalization 
cannot occur over such short times (AT thcrmal ~ 0.5 s), and since the neutrino production rate is sensitive to the local 
temperature at the production site, the inhomogeneities during the collapse phase get reflected in the inhomogencities 
in the escaping neutrino fluxes and their distributions. Because of both the v spin-magnetic field (B) and v spin- 
angular momentum (J) coupling the asymmetries in these distributions can give rise to quadrupole moments, which 



1 The attentive reader must regard that neutrinos carry away almost all of the binding energy of the just-born neutron star, 
i.e., A£„ - 3 x 10 53 erg. 
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must generate gravitational waves as suggested by Mosquera Cuesta (2000; 2002), and dipole moments which can 
explain the origin of pulsar kicks (Kusenko & Segre 1997; Fuller et al. 2003). 

Fixed by the probability of oscillation, P^Z^^ , the fraction of neutrinos that can escape in the first few milliseconds 
is, however, small. Firstly, the neutrinos have to be produced roughly within one mean free path from their resonance 
surface. Secondly, since in the case of v a — ► v s oscillations m s is the heaviest neutrino species, the sign of the effective 
potential A(x) (see discussion below) and the resonance condition indicates that only ^ e s and the antineutrinos 
and v T can undergo resonant conversions. In Section IV we address all these issues and determine this fundamental 
property of the mechanism for producing GWs from v flavor conversions. 

III. ANISOTROPIC NEUTRINO OUTFLOW: ORIGIN AND COMPUTATION 

To provide a physical foundation for the procedure introduced here to determine the neutrino asymmetry parameter, 
a, which measures how large is the deviation of the v flux from a spherical one, we recall next two fundamental effects 
that run into action once a PNS is forming after the supernova collapse. We stress that other physical process such 
as convection, thcrmalization, etc., are not relevant over the time scale under consideration: AT osc ^3—10 ms after 
the SN core bounce. Those effects take a more longer time (<~ 30 — 100 ms) to start to dominate the physics of the 
PNS, and therefore do not modify in a sensitive manner the picture described below 2 . Indeed, if v thermalization, 
for instance, already took place, then the oscillations are severely precluded since oscillations benefit of the existence 



2 At this point, a note of warning regarding the time scale we are using for the present calculations is of worth. This is specially 
so in the light of the very recent paper by Loveridge [24], where a very important extension of our original idea, introduced in 
Refs. [28], is detailedly provided. In his computation of the gravitational radiation emission from an off-centered flavor-changing 
v beam Loveridge used a time scale of /\j ll/ - burst ^ iq Sj apparently based on the duration of the i^-burst from SN1987A. 

The great novelty in Loveridge's paper is the prediction of a periodic GWs signal from the flavor-changing v beam eccentrically 
outflowing from the just-born pulsar. This GWs signal would have characteristics that make it observable by both LIGO and 
LISA GWs interferometers. Interestingly in itself, various arguments favouring evidence for regularly pulsed neutrino emission 
from SN1987a in the period range between (8.9-11.2) s were given by Harwit, et al. [18] and by Saha and Chattopadhyay [35]. 
However, Fischer checked systematically for periodicities between (5-15) ms in the v burst from SN1987A [13]. He disclaimed 
all those hypotesis by showing that the multitude of medriocre period fits seemed to be rather typical for events distributed 
randomly instead of periodically [13]. Thus no significant periodicity exists in the arrival time of the neutrinos from SN1987A 
as detected by Kamiokande II and 1MB detectors. Nonetheless, any evidence for such a regular v signal in a forthcoming 
(future) supernova might decidedly favor Loveridge's GWs mechanism from v oscillations in nascent pulsars. 

As it stands, however, Loveridge's mechanism is decidedly different from ours in several respects. Firstly, we do not invoke an 
off-centered rotating v beam for producing the GWs emission from the v oscillations. In our mechanism the neutrino outflow 
is simultaneously acted upon by both the pulsar centered magnetic field and angular momentum vectors, as we describe below. 
Thus, the ^-spin coupling to both vectors turns out to be the source of the, at least, quadrupolar v outflow and GWs emission 
during the flavor conversions. Secondly, as a consequence of this v escape from the star the resulting GWs signal, in principle, 
is not periodic, as opposed to Loveridge's. Indeed, the GWs signal will look much like the one computed by Burrows and 
Hayes [7], including the appearance of a Christodoulou's memory in the waveform. Thirdly, as discussed next, the overall time 
scale for the process to take place in our mechanism is about three orders of magnitude shorter compared to one assumed by 
Loveridge. 

The v oscillation time scale and the duration of the GWs emission are both crucial features of both the mechanisms above 
discussed. In this regards, we advise that a very extensive set of references (here we just quote a few of them) showed that the 
v signal from SN1987A exhibits a peculiar time profile [37,3,9]. According to Refs. [37,3,9], the v burst observed from SN1987A 
is bunched into three clusters around (0.0-0.107) s, (1.541-1.728) s, (9.219-12.349) s. In particular, Cowsik [9] claimed that 
during the very early phase the lighter of the neutrinos arrived and by AT ~ 0.1 s all of the neutrinos above the 1MB threshold 
of 20 MeV had already gone past and thus were not seen by this v detector. Therefore, if one stands on these pieces of evidence 
it is clear that the largest portion of the total ^s from SN1987A were emitted over a time scale smaller than 100 ms. This last 
time scale is an order of magnitude larger than the one we are favouring here, one which is taken from most of the theoretical 
analysis and numerical simulations of supernova core collapse and v emission (see for instance Refs. 10,11,12 in Sato and Suzuki 
[37]; and also Burrows and Lattimer [6]; and Refs. [43,25]). A possible explanation of such a difference may be that a large 
portion of the released electron v& undergone flavor conversions into sterile vs over a time scale much shorter than 100 ms, a 
reason of why they went undetected. 

In brief, although from the observational point of view, we may agree that the (~ 10 s) v emission time scale assumed by 
Loveridge [24] appears a reasonable one, we think of it is not a very realistic one to picture out the v flavor conversion mechanism 
inside the nascent neutron star, specially if one takes into account that the v oscillation process implies a very short time scale: 
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of energy, matter density and entropy gradients inside the PNS (Bilcnky, Giunti and Grimus 1999; Akhmedov 1999), 
which are "washed out" once thermalization onsets. 

EDITOR PLEASE PLACE TABLE 1 AROUND HERE !!! 



A. Why no room for ^-driven convection over AT os 



In order to back the dismissal in our discussion on neutrino oscillations of the effects of convection inside the proto- 
neutron star, we would like to take advantage of some arguments presented in the state-of-the-art of the subject by 
Janka, Kifonidis & Rampp (JKR, 2001), who provided detailed analysis about convection inside the nascent neutron 
star. In particular, these authors showed that the growth time scale of convective instabilities (r cv ) in the neutrino- 
heated region (adjacent layers outside the just-born neutron star, of relevance for successful supernova explosions) 
depends on the gradients of entropy and lepton number through the growth rate of Ledoux convection, ox, as 
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or equivalently 
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where Mi = 1M Q , and R s and R g define the shock and gain radius, respectively, and R g j is a function of the gain 
radius, the temperature inside the star, and the neutrino luminosity and energy. Here the estimates were obtained 
for g — GM/Rg, (dp/ds) P ~ —p/s, and ds/dr ~ —\ ^ R S _ R ^ (see JKR 2001, for further details). 

Numerical simulations demonstrate that convection inside the proto-neutron star does start as early as a few tens 
of milliseconds after core bounce. It develops in both (i) unstable surface-near regions, i.e., in layers around the 



neutrinosphcre where density is p < 10 g cm" 



and (ii) deeper layers of density p > 10 12 g cm 3 ) 



where a negative 



lepton number gradient appears. Despite of this piece of evidence, the time scale defined by Eq.(5) is relatively long 
compared to both the estimated time interval for the deleptonization process to take place: ATdoiopt ~ 10 ms, i.e., 
time over which most electron neutrinos are produced (Mayle, Wilson & Schramm 1987; Walker & Schramm 1987), 
and the core bounce time scale: ATbouncc ~ 20 ms; where the large part of the neutrino luminosity associated with 
other flavors is produced through processes like bremsstrahlung, neutrino-neutrino and neutrino-nucleon scattering 
within less than 5 ms (Mayle, Wilson & Schramm 1987; Walker & Schramm 1987; JKR 2001, and references therein). 
Indeed, from the convective regions below the neutrinosphere neutron fingers dig into the star and reach its center 



the resonance or coherence length time scale. As claimed above, this time scale must be related to the oscillation length (not the 
system's response time scale) over which most of the conversions must take place. Moreover, oscillations of massive neutrinos 
are damped when the propagation distance is greater than the coherence length 

icoh = A osc y'^j^j ■ ( 3 ) 

Here E and AE are, respectively, the v energy and energy spread determined by the production and detection conditions. In 
a supernova, the neutrinos nonforward scatter in continuous energy distributions so that AE ~ E, and hence the coherence 
length is nearly the oscillation length [32], which fixes the time scale AT oac we call for henceforth. Yet, in the early phases of a 
SN the neutrino flux is so large that the weak-interaction potential created by the neutrinos is compatible to that of the baryon 
matter around. Thus, neutrinos can be thought of as a dominant background medium that acts as a coherent superposition 
of flavor states that drives the conversions in a nonlinear way. In other words, the oscillations become "synchronized" , which 
means that all modes oscillate collectively with the same frequency [36,33]. Such a frequency should be related to the oscillation 
or coherence length, and through it to the oscillation time scale. Thus, this last behaviour adds to our argument in favour of 
a shorter time scale for the overall oscillations to take place. So, the AT osc < 10 ms is well-fundamented. Besides, a time scale 
that long as Loveridge's ~ 10 s strongly disagree with Spruit and Phinney constraint on the overall time scale: ~ 0.32 s, for 
the kick driving mechanism [38]. 
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in about one second. Then they propagate outwards to englobe almost all the exploding star. Under the physical 
conditions dominant over the first 10-20 ms after core bounce one expect most neutrino oscillations of all flavors to 
take place at that time. Thence, convective effects are not relevant during a time scale that short. As such, it cannot 
modify in a significant fashion our analysis regarding the mechanism for the generation of gravitational waves from 
neutrino oscillations, here highlighted. 

Moreover, JKR (2001) stressed the "desastrous" role of rotation for convection. A high rotation velocity of the just- 
born neutron star reduces dramatically the effects of convection because of the suppression of the neutrino-nucleon 
interaction due to nucleon correlations in the nuclear medium composing the proto-neutron star. Physically, rotation 
leads to a suppression of convective motions near the rotation axis because of a stabilizing stratification of the star 
matter specific angular momentum. In passing, we stress that a similar effect is also expected from the action of a 
background magnetic field. Both effects, rotation and magnetic field, then appear to be more crucial for the physics 
of neutrino interactions inside the newly-born neutron star, and for the production of GWs during the oscillation 
transient. Thus we address both of them next. 

B. ^-rotation interaction 

That gravity couples to neutrinos is well-known since Dirac. The very first work, as far as we are aware of, to 
show that the particle spin and PNS rotation couples gravitationally, a universal feature (Soares & Tiomno 1996), 
was that of Unruh (1973). It was shown that a consistent minimal-coupling generalization to a curved background 
(a Kerr spacetime in that case) of the v equations is possible, and that it leads to equations separable for the radial 
and angular components, though coupled. For a massless v field (compare to Eq.(2)), and standard Minkowski space 
Dirac matrices r y A , the Dirac equation derived by Unruh reads 

7^-1^ = 0, (6) 
where the Dirac 7^ matrices relate to the Kerr spacetime metric through 

1 A 1 B +1 B 1 A = 2g AB^ (?) 

Eq.(7) shows that the j A matrices satisfy the Clifford's algebra. Further, the spin-affine connections Ta in Eq.(6) 
are uniquely determined by the relations 

T A1 B - 7 B r A = -JL + T aA V , and ti(T A ) = . (8) 
ox A 

From the ^-number current, J (x) — ^(x)j tp(x), it was shown that the {/-number density is always positive and 
given as (—g) 1 / 2 J t (x) = (—g) 1 / 2, ip(x)~/ t ip(x) (see Unruh 1973 for details). The complete analysis of the coupling shows 
that the v field in this background is not superradiant, as opposed to the case of the classical fields studied previously. 

Vilenkin (1978) extended the above analysis and showed that, upon admitting helicity (L) to be a good quantum 
number, the angular distribution Fj m (9) of the thermal fermion gas of ^s (L = +1) and ^s (L = —1) in the mode (j, m), 
with the function Fj m (9) satisfying the normalization condition 2tt Fj m (9) sin 6 d6 = 1, leads to an asymmetric v 
emission (see illustration in Fig. 2) from a Kerr black hole (BH), of specific angular momentum a = J/M, described 
by 



dN 1 22 J2±( l ± Lac ° s0 ) 

dtdudO ~ 8tt2 ^ f e [f>±ifiBH)] + 1 



= ^MW T ^— . '- . (9) 



Eq.(9) shows that more us are emitted in the direction parallel to the BH's spin, whilst more vs escape in the 
antiparallel direction. 3 Further, for other weak-interacting particles emitted from the BH parity is not conserved 
(Vilenkin 1978). More fundamental yet, Vilenkin (1978) work demonstrates that the same physics must be valid for 



3 Notice that this behaviour is also manifest in the case of the magnetic field i^-spin coupling discussed below. Therefore, these 
effects affect the oscillation probability, as we discuss later on. 
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any other rotating star. In other words, the v spin coupling to rotation, in a gravitational background, is a universal 
feature, regardless which the spacetime source can be. Here onwards we shall take advantage of this feature for the 
case of neutrino emission from a rotating NS, and suggest that the basic quadrupole nature of the v emission from 
the PNS stems in part from this spacetime effect. The other fundamental effect we address right next. 

EDITOR PLEASE PLACE FIGURE 2 AROUND HERE !!! 



C. v-B field interaction 



That the electromagnetic properties of vs are modified due to its interaction with a background matter distribution 
is a well-known fact. This is reflected in the additional contribution to its self-energy stemming from the v spin-to- 
magnctic field coupling (B ■ k v ). In the case of v T — > v e oscillations, for instance, the geometry of the ^-sphere is 
dramatically deformed by B (see Fig. 2), an effect that strongly depends upon the relative pointing directions of both 
B and v momentum k v . The magnetic coupling distorts the ^-surface ("sphere") in such a way that it is no more 
concentric with the ^ e -sphere (see Fig.l and 2 in Kusenko 1999). Therefore, v T s escaping paralelly-pointing to the B 
field have a lower temperature than those flowing away in the opposite (anti-parallel) direction. 

At the same time, the v spin coupling to rotation, Eq.(9) (Vilcnkin 1978), drives also an effective momentum 
(and thus energy flux) asymmetry along the angular momentum direction J, as shown in Fig. 2. For a relative 
orientation 9(B «-> J) ^ between B and J, i. e., a canonical pulsar, this combined action on the escaping ^s of a 
rotating background spacetime plus magnetic field makes their i/-sphere a decidedly distorted surface. More precisely, 
the volumetric region obtained by rotating around the hatched regions in Fig. 2, becomes at least a quadrupolar 
outflowing energy distribution. This is the source of the strong GWs bursts in this mechanism when the oscillations 
ensue. It also justifies the large value of the anisotropy parameter a used in Eq.(14) below (see also Burrows & Hayes 
1996). 

EDITOR PLEASE PLACE TABLE 2 AROUND HERE !!! 



D. The anisotropy parameter 

Based on the concommitant action of both effects: the v-sp'm coupling to both the magnetic field and rotation 
described previously, one can determine the v flow anisotropy in a novel, self-consistent fashion by defining a as the 
ratio of the total volume filled by the distorted ^-spheres to that of the proto-neutron star (PNS), as one can infer 
from Fig. 2. The ^-sphere radius of a non-magnetic non-rotating star is obtained from the condition: 

T Ve (R»e) = j R £»e P(r) dr =l> ( 10 ) 

where T„ e is the optical depth, and K Ve the scattering opacity for electron neutrinos, and p the matter density. 
Following Burrows, Hayes & Fryxell (1995) one can take hereafter R Vc = i?p e NS = R v ~ 35 km, which is of the order 
of magnitude of the oscillation length A„ of a typical supernova v, constituent as well of the atmospheric ^s for which 
Am 2 ~ 10 _3 eV 2 has been estimated by Superkamiokande v detector (Fukuda et al. 1998) 



A„ ~ 31 km 



E u 



10 MeV 



in- 3 eV 2N 

— — ) (ID 

Am 2 ' { ' 



Therefore, resonant conversions between active species may take place at the position r from the center defined by 

r = R Ve + So cos <f) , (12) 

with 6 the angle between the ^-spin and B, i.e., cos0 = ^ k '^ , and 

k 

*-5?fc~ 1 - Mkm - (13) 
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for B ~ 10 14 15 G, respectively. Here e, N e = Y e N n and fi e represent the electron charge, density and chemical 
potential, respectively. This defines in Fig. 2 an ellipsoidal figure of equilibrium with semi axes 

a = R v + So; and b = R v , (14) 

and volume (after rotating around B): 

V ellips . = ^kRI{R v + 8 ) . (15) 

Meanwhile, the i^-spin coupling to rotation described by Eq.(9) generates an asymmetric lemniscate-like plane curve 
(see Fig.2) 

r = R v (1 ± L a cos 9) , (16) 
which upon a 2tt rotation around the star angular momentum axis generates a volume: 

V lemnisc . = - A {2R v f x 2ny + \{Ruf x 2iry , (17) 



where the quantity y (and x) is defined as the location of the centroid of one of the lobes of that plane figure with 
:spect to its coordenate center (x,y), a 
straightforward, calculation one obtains 



respect to its coordenate center (x,y), and is obtained from the standard definition y — j^jj- ■ After a long, but 



v = f^i) Rv ' *=^m R - (18) 

Thence, for a PNS as the one modeled by Burrows, Hayes & Frixell (1995), with parameters as given in Table II 
one obtains 

Vellips. Vlemnisc. r\ t t r\ r\t /-ir\\ 

a min = £— 0.11 - 0.01 , (19) 

^PNS 

a figure clearly compatible with that one in Burrows & Hayes (1996). As such, this is essentially a new result of 
this paper. The attentive reader must notice in passing that the definition in Eq.(19) does take into account all of 
the physics of the neutrino oscillations: luminosity, density gradients and angular propagation, since Eqs.(12,16) do 
gather the relevant information regarding the spatial configuration of the v luminosity in as much as is done in the 
standard definition of a (Burrows & Hayes 1996; Miiller & Janka 1997) 

< *>-z^i*' , <™ i # fl - (20) 

Indeed, one can get the "flavour" of the relationship between these two definitions by noticing that the quantity 
L„(n',t) in the integrand of Eq.(20) can be expressed as L V (Q! ,t) = L v (t)F(Cl'), where the function F(Q') contains 
now all the information regarding the angular distribution of the neutrino emission. Hence L v {t) can be factorized out 
of the integral and dropped from Eq.(20). This converts Eq.(20) in a relationship among (solid) angular quantities, 
which clearly can be reduced to a volumetric one, similar to the one introduced in Eq.(19), upon a transformation 
using the definition of solid angle in the form of Lambert's law. e?f2 = £p s , and applied to the sphere representing 
the PNS. Here 9, measured from a coordinate system centered at the PNS (source coordinate system in Miiller & 
Janka 1997), plays the role of the angle between the direction towards the observer and the direction Q' of the 
radiation emission in Eq.(24), and subsequents, in Miiller & Janka (1997). Therefore, the novel result here presented 
is physically consistent with the standard one for the asymmetry parameter a. 



IV. ENLARGED v AND GWS LUMINOSITY FROM OSCILLATIONS IN DENSE MATTER 

v outflow from a SN core bounce is a well-known source of GWs ( Epstein 1978; Burrows & Hayes 1996; Miiller 
& Janka 1997; Mosquera Cuesta 2000; 2002). Numerical simulations (Miiller & Janka 1997) showed that, in general, 
the fraction of the total binding energy emitted as GWs by pure v convection is: E^ w ~ [10~ 10 -10~ 13 ] M Q c 2 , for a 
v total luminosity: L v ~ 10 53 erg s _1 . 
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Unlike GWs produced by v convection (Miiller & Janka 1997), in the production of GWs via v oscillations (Mosquera 
Cuesta 2000; 2002) (y a <-> v a or v a «-> v s ) there exists two main reasons for expecting a major enhancement in the 
GWs luminosity during the transition: a) the conversion itself, which makes the overall luminosity (L Vx ) of a given v x 
species to grow by a large factor: AL Ux < 10%-L^ ^, see below. The enhancement stems from the mass-energy being 
given to, or drained from, the new v species into which oscillations take place. This augment gets reflected in the 
species mass-squared difference, Am 2 , and their relative abundances: one species is number-depleted while the other 
gets its number enhanced. But, even if the energy increase, or give up, is small, b) the abrupt resonant conversion 
over the transition time (see Table I) 

ATosc = T -/° SC ~ [10- 4 - 10- 3 ] s, (21) 

K-Diff. 

also magnifies transiently L„ x . Here Aosc defines the oscillation length (computed below), and V v -d% ff. ~ 10 9 cms -1 
the convective v diffusion velocity (Miiller & Janka 1997). In Section IV we estimate the transition probability: 
P Ua ^ Ua (\x — xq\), P Va ^ Vs (\x — xq\), the quantity that measures how many vs can indeed oscillate. This probability 
also fixes the total amount of energy participating in the generation of the GWs through this mechanism, as shown 
in Section V. 

If flavor transitions can indeed take place during supernovae (SNe) core collapse and bounce, then they must leave 
some imprints in the SNe neutrino spectrum. Main observational consequences of neutrino conversions inside SNe 
include a) the partial or total disappearance of the neutronization peak; the moment at which most v e s are produced, 
b) the interchange of the original spectrum and the appearance of a hard v e spectrum, together with c) distorsions of 
the v e energy spectrum and d) alterations of the v e spectrum (Dighe & Smirnov 2000) . As discussed below in Section 
4.4, observations of the neutrino burst from SN1987a have allowed to put some bounds on both v a — > v a and v a — > v s 
classes of flavor conversions. 



A. Resonance, adiabaticity and the role of weak magnetism 

As argued by Mosquera Cuesta (2000; 2002) v oscillations in vacuum produce no GWs. In the case of active-to- 
active v oscillations (essentially the same physical argument holds also for active-to-sterile v oscillations), the main 
reason for this negative result is that this class of conversions do not increase in a significant figure the total number 
of particles escaping from the proto-neutron star. In the case of active-to-activc v oscillations in dense matter, the 
process generates no GWs since the oscillations develop with the neutrinos having very short mean free paths, so that 
they motion outwards can be envisioned as a standard difussion process. 

However, if one takes into account the novel result by Horowitz (2002) the situation may change dramatically. 
According to this author, because of the active antineutrino species weak magnetism their effective luminosity can be 
enlarged as much as 15% compared to the typical one they achieve when this effect is not taken into consideration 
during their propagation in dense matter. This result can be interpreted by stating that the number of oscillating 
(and potentially escaping) antineutrinos may be augmented by a large factor because now they do encounter longer 
mean free paths. Below we take advantage of this peculiar behaviour of v outflow in supernovae in computing the 
overall probability of transition between active species and the GWs emitted in the process. 

EDITOR PLEASE PLACE TABLE 3 AROUND HERE !!! 

Wolfenstein (1978; 1979) and Mikheev & Smirnov (1985) pointed out that the neutrino oscillation pattern in 
vacuum can get noticeably modified by the passage of neutrinos through matter because of the effect of coherent 
forward scattering. Therefore, interaction with matter, as pictured by Eq.(l), may help in allowing more ^s to escape 
if resonant conversions into active (Walker & Schramm 1987) and/or sterile vs (Mosquera Cuesta 2000) occur inside 
the ^-sphere of the active vs. 

The description of the two-neutrino oscillations process in matter follows from the Schrodinger-like (because the 
dynamics is described as a function of the space variable x instead of the standard time t) differential equation 
(Bilenky, Giunti & Grimus 1999; Grimus 2003) 

. d / a a \ = t / a a \ 
dx V a P J V a & J 
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AE 



i\+ml + 2^G F {N{ Va ) + N(v p )) 







i 





+ 



where 



A - Am 2 cos 29 Am 2 cos 29 
Am 2 cos 29 -A + Am 2 cos 20 



A = 2\/2GfE (N(u a ) — N{vp)) , 



:)}• 



(23) 



and N(v a ) = 5 ae N e — \N n , a, (3 — e,/j,,T,s, N(v s ) = 0, and Am 2 = m\ — m\. The eigenfunctions of the matter 
effective Hamiltonian follow from the relation 



where 



ipmi = 



cos f, 
— sin( 



1pm2 = 



sin t 
cos( 



(24) 



(25) 



The eigenvalues of Ej and the matter mixing angle 9 m are thus given as 



-El,: 



m 2 + m 2 + 2V2G F (N(is a ) + N{vp)) 
T yf{A - Am 2 cos 29) 2 + {Am 2 sin 29) 2 



and 



tan2#„, = 



tan26» 



1 - 



(26) 



(27) 



Am 2 cos IB 



where 9 is the vacuum mixing angle. By defining U m {x) = (ip m i, ipmn) as the mixing matrix of n v, and 
Sj = —i f* 1 dx' tpmj{x'y^p m j(x') as the adiabatic phases, one can compute the oscillation amplitude as 



5, + 



dx'Ej(x') 



X U^{x )c 



(28) 



Finally, by averaging over neutrino energies, i.e., by setting (exp {—i dx' AE(x')j'^ = (sec Bilenky, Giunti 
& Grimus 1999), the outcoming transition probability among active flavor-changing species thus reads 



p va^va(\xi - ^o|) = 1 - 2 I 1 + cos26' m (a;o)cos26' m (a;i)] , 
where x and X\ correspond to the production (emission) and detection sites, respectively. 



(29) 



B. Active-to-active v oscillations 



As stated above, in order to produce an effect neutrinos must be able to escape the core without thermalizing with 
the stellar material. For active neutrino species of energies fO MeV, this is not possible as long as the matter 
density is > fO 10 g cm~ 3 . Since the production rate of neutrinos is a steeply increasing function of matter density 
(production rate oc p n , where p is the matter density and n > f), the overwhelming majority of the neutrinos of 
all species produced are trapped. This way, there seems to be no contribution to the GWs amplitude for neutrino 
conversions taking place within the active neutrino flavors. In the first paper of this series (Mosquera Cuesta 2000), 
this difficulty was overcome only by addressing neutrino conversions into sterile species. Nonetheless, if there were 
indeed weak magnetism effects (Horowitz 2002), one can rethink conversions within active species. In his new result on 
weak magnetism for antineutrinos in core collapse supernovae, Horowitz (2002) showed that the antineutrinos (v x s) 
luminosity could be noticeably increased because of their longer mean free paths, and this means that the total energy 
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flux can be augmented in ~ 15% for us of temperature ~ 10 MeV. One can verify that longer mean free paths allows 
for a larger oscillation probability, and hence the contribution to the generation of GWs during flavor conversions 
within active species becomes nonnegligible compared to the earlier case (Mosquera Cuesta 2000) where not weak 
magnetism effects were taken into account. 

For active-to-active v flavor conversions, for instance: v e — > u^, v T \ as implied by SNO results, the resonance must 
take place at a distance x YCS from the PNS center and amid the active ^-spheres, whenever the following relation is 
satisfied (see Table III) 



Am 2 cos2(9 m = 2\/2 G F N e (x res ) k Ve , 



(30) 



notice that we neglected the magnetic field contribution. Here k Ve — E Uc /c is the v e momentum, and N e (x Tes ) = 



10 39 40 cm 3 is the electron number density. Thus the right-hand part of this equation reduces to 



2V2G F E Ue N e (x ies ) = 6.9 x 10 4 eV 2 



10 11 g cm 3 



10 MeV 



(31) 



For densities of order p x 



10 u g cm 3 , i.e., recalling that v e s arc produced at the PNS outermost regions 



(Walker & Schramm 1987) where the electron to baryon ratio is Y e ~ 0.1, the resonance condition if satisfied for a 
mass-squared difference of about Am 2 ~ 10 4 eV 2 , which implies a neutrino mass of about ra v <~ 10 2 eV. Neutrino 
flavor conversions in the resonance region can be strong if the adiabaticity condition is fulfilled (Walker & Schramm 
1987), i.e., whenever (Bilenky, Giunti & Grimus 1999; Grimus 2003) 



Am 2 sin 2 20 n 
2E V cos 29 m 



1 dp 

p dx 



» 1 



(32) 



where x res is the position of the resonance layer. Recalling that the typical scale of density variations in the PNS 
core is /ijv e ~ (dNg/dr) -1 ~ 6 km, this adiabatic behaviour could be achieved as far as the density and magnetic field 
remain constant over the oscillation length 



-^osc — 



/ 1 dp 



/ 1 Am 2 



P J x=Xres \2ir 2k v 



sin(26» m ) 



1 cm 
sin(20 TO ) 



of order hjy e , which can be satisfied for Am 2 <~ 10 4 eV 2 as long as (Kusenko & Scgre 1999) 

sin 2 29 m > 1Q- 8 . 



(33) 



(34) 



Although these v oscillations could be adiabatic for a wide range of mixing angles and thus a large number of 
^s could actually oscillate, a v x mass such as this is incompatible with both viable solutions to the Solar Neutrino 
Problem (SNP) and the most recent cosmological constraints on the total mass of all stable neutrino species that could 
have left their imprint in the Cosmic Microwave Background Radiation (CMBR), as inferred from the observations 
performed by the satellite WMAP: m v <~ 1 eV. Therefore, we dismiss this possibility since there appears to be no 
evidence for neutrinos masses in this parameter range (m Vx <~ 10 2 eV) inside the PNS core. 

On the other hand, if one takes into account the KamLAND results (Eguchi et al. 2003), which definitively 
demonstrated that 

i) a large mixing angle (LMA) solution of the solar v problem is favoured: sin 2 26 ~ 0.8, 

ii) for a mass-squared difference: Am 2 <~ 5.5 x 10~ 5 eV 2 (we use the approximate value Am 2 ~ 10~ 4 eV 2 for the 
estimates below), 

one can see that resonant conversions with Am 2 = m 2 2 — <~ 10 -4 eV 2 would take place in supernova regions 
where the density is about p <~ 10 — 30 g cm~ 3 , which corresponds to the outtermost layers of the exploding star. 
Although a large number of v species can in effect participate of the transitions there, i.e, the neutrino luminosity 
can be still a very large quantity, these regions are of no interest for the gravitational-wave emission from neutrino 
oscillations since the overall energy density at that distance from the star center is relatively small. This does not 
mean that no GWs are emitted from transitions there, it is to mean that their strain is very small so as to be 
detectable. As discussed by Dighe & Smirnov (2000), observations of v oscillations in this parameter range would 



provide useful information regarding the SNP, the hierarchy of neutrino masses and the mixing | U, 
that oscillations in this range would imply a mass for the vi species m„ 2 <~ 10~ 2 eV 



■ 10 g cm" 



e 3| _ . Note in passing 
r), in the case when 



> m Vl . This m„ 2 is compatible with current limits from WMAP (Pierce & Murayama 2003; Hannestad 2003). 
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Finally, let us consider oscillations in the parameter range estimated from CMBR by WMAP observations. In this 
case, resonant v transitions would take place in regions where the density is as high as p ~ 10 7 ~ 8 g cm ~ 3 , that is, at 
the supernova mantle or PNS upper layers. At these densities the oscillation length can be still A osc ^1 — 5 km, and 
thus the conversions can be considered as adiabatic. Thus the resonance condition can be satisfied for Am 2 ~ 1 eV 2 
as long as sin 2 29 m < 10~ 3 . 

Hence, by plugging this constraint into Eq.(27), and recalling that i) at least 6 ^-species can participate in the 
flavor transitions, ii) most vs are emitted having fc„ parallel to B, which implies a further reduction factor of 2, and 
also Hi) most ^s are emitted having k v parallel to J implying an additional reduction factor of 2, one can show from 
Eq. (29) that the fraction of v a species that can eventually exchange flavor during the first few milliseconds after core 
bounce turns out to be (Mosquera Cuesta 2000) 

i , „ a ^„ a (|j?-2rea|)~0.1%. (35) 

Thus, the total energy involved in the oscillation process we are considering could be estimated as: £7* ot ~ 
P Ua ^ Ua (\x — x res \) x N v k v c (1 + 0.15), with N v the total number of neutrinos undergoing flavor conversions during 
the time scale AT osc . This value leads to a bit stronger GW burst, as we show in Section V below. 

Note in passing that KamLAND v e experiment suggests fp e ->p ~ 40%! (Eguchi, et al. 2003), while LSND 
Pv^v c ~ 0.26%. This result from LSND has not been so far ruled out by any terrestrial experiment, and there is a 
large expectation that it could be verified by MiniBoone at Fermilab. 

EDITOR PLEASE PLACE TABLE 4 AROUND HERE !!! 



C. Active-to-sterile v oscillations 



On the other hand, if the Liquid Scintillator Neutrino Detector (LSND) were a true indication of oscillations 
(P M <-» v e ) (Pierce & Murayama 2003; Hannestad 2003) with Am 2 SND <~ 10 _1 cV 2 , then active-to-stcrile v a — > v s 
oscillations could take place during a supernova core bounce (see Mosquera Cuesta 2002). Such oscillations would be 
resonant whenever the resonance condition is satisfied. As seeing from Table III this could happen for 



V2G F 



N e (x)--N° n (x) 



A 777.^ 

= V(x) = -—cos 2^ 
ihi,. 



(36) 



Table III also shows that for f MiT <— 



v s the N e term is absent, while in the case of vs, the potential V(x) changes 
by an overall sign. Numerically, for v e <-> v s oscillations 



V(x) 



7.5xl0 2 ( cV 



MeV 



10 10 g cm" 3 



SYe 

2 



(37) 



For v^ T <-> v. 



oscillations, the last term in parenthesis becomes (-rf — h) ~ 1 



v conversions in the resonance region prove to be enhanced if the adiabaticity condition is fulfilled (Walker & 



Schramm 1987). This is the same as requiring the oscillation probability in Eq.(29) to become P„ a 



cos 2 29 Tl 



1. 



Moreover, after the resonance region the newly created sterile ^s have very a small probability (< Pv e -w a >= 
\ sin 2 29 m ) of oscillating back to active ^s, which could be potentially trapped. It is easy to check that the resonance 
condition in Eq.(36) is satisfied whenever 



1 eV 2 < Am 2 < 10 4 eV 2 



(38) 



Meanwhile, the adiabaticity condition, Eq.(32), holds if 

, sin 2 29. 



Am 2 



cos 29, n 



> 10" 3 eV 2 



lOMeV 



(39) 



since at the PNS core A osc ~ 1 km (see Eq.(33)). This is easily satisfied for Am 2 < 10 4 eV 2 as long as sin 2 29 m < 
10- 7 . This limit on the mixing angle is in agreement with Superkamiokande strong constraints on v e — > v s mixing 
(Nunokawa 2001). Thence, we find that a substantial fraction (P„ a ^„ s ~ 1%) of ^s may get converted to sterile vs, 
and escape the core of the star, if the sterile v mass (m U4 ) is such that 
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1 eV 2 < Am^„ 4 < 10 4 eV 2 



(40) 



By looking at the lower limit for Am 2 in Eq.(40) one can say that such a mass difference cannot solve the observed 
solar v problem and/or being compatible with the atmospheric v observations, but the possibility of having three light 
active vs of mass m Vl 2 3 <~ [10~ 2 — 10~ 3 ] eV explaining these anomalies and a "heavy" sterile v s of mass m„ 4 < 1 eV as 
required by the LSND experiment remains to be a viable alternative. On the other hand, if the sterile neutrino could 
be as massive as m Vi < 1 keV this mass will make it a very promissing candidate as a constituent of the universe's dark 
matter. This last possibility was recently readdresscd by Fuller et al. (2003), who estimated sterile neutrino masses 
in the range (1-20) keV and small mixing angle sin 2 29 <~ 10 -4 with the electron neutrino as a potential explanation 
of pulsar kicks. The same argument had also been considered by Mosquera Cuesta (2002); and references therein. 

As quoted above, for both classes of v conversions the number of ^s escaping and their angular distribution is 
sensitive to the instantaneous distribution of production sites. These inhomogeneities can give rise to quadrupole 
moments that generates GWs (Mosquera Cuesta 2000; Miiller & Janka 1997), and dipolc moments that could drive 
the runaway pulsar kicks (Kusenko & Segre 1997; Fuller et al. 2003). Noting that at least 6 ^-species can participate 
in both types of oscillations and that the interaction with both the magnetic field and angular momentum of the PNS 
brings with an overall reduction factor of 4 in the oscillation probability, one can show that the fraction of vs that 
can actually undergo flavor transitions in the first few milliseconds is (Mosquera Cuesta 2000) 

p^:(\x-x Q \)<m, (4i) 

of the total ^s number: M Q x 1 ~ 10 57 V, which corresponds to a total energy exchanged during the transition 

AE^2Xl s ~3x 10 51 erg. (42) 

Equations (42,21) determine the total luminosity of the neutrinos participating in the resonant flavor conversions. 
These figures are called for in the definition used in Eq.(43) below, as the basis to compute the characteristic amplitude 
of the GWs emitted during the transitions discussed above. 

The attentive reader must notice, however, that the lower limit in Eq.(38) stems from using the constraint derived 
from the observations of the CMBR performed by the satellite WMAP, which suggested that the total mass of all 
neutrino species should not be larger than 1 eV. In a hierarchy where the heaviest v is the sterile v s , this limit leads 
to a maximum mass: m™ ax <~ 1 eV. The use of this mass-squared difference constraint implies that the density at the 
PNS region where the oscillations can take place is p^ MAP ~ 10 8 g cm" 3 . If the WMAP constraint on m v stands up 
also for sterile neutrinos, this would preclude the mechanism for giving to nascent pulsars the natal velocities (kicks) 
from active neutrino flavor conversions into sterile neutrinos, as claimed by Kusenko & Segre (1997) and Fuller et 
al. (2003), since in such a case the transitions would take place in regions far outside the PNS core from where no 
influence could be received back once the oscillations develop. Nonetheless, the GWs emission would still take place, 
as we stressed above. 

D. Experimental bounds on supernova v oscillations 

Studies of SN physics have also focused on the potential role of oscillations (Walker & Schramm 1987) between active 
and sterile vs. In particular, there are limits on the v e <-» v s conversion rate inside the SN core from the detected v e 
flux from SN1987a (Nunokawa et al. 1997; Nunokawa 2001). According to (Nunokawa et al. 1997; Nunokawa 2001), 
the time spread and the number of detected v e events constrain v e <-> v s oscillations with: 10 6 < Am 2 < 10 8 cV 2 for 
10 -3 < sin 2 20„ e _ >1/> < 10~ 7 . More stringent constraints stem from arguing that if there were too many "escaping 
vs" , the SN explosion itself would not take place (Nunokawa et al. 1997; Nunokawa 2001). Such bounds are, however, 
model dependent. One should keep in mind that the mechanism through which the explosion takes place is, in fact, 
not well established (Buras, Rampp, Janka & Kifonidis 2003). In effect, it has been recently claimed that even after 
including all the best physics we know today about, i.e., the state-of-the-art on supernova physics, the numerical 
models of exploding stars do not explode as they should (Buras, Rampp, Janka & Kifonidis 2003). There seems to 
be a piece of missing physics in those formulations. Thence, there seems to be no hope of achieving P^ a ~^ a larger 
than <~ 1% during v oscillations in supernovae. 
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V. GWS ENERGETICS FROM v LUMINOSITY AND DETECTABILITY 



If v oscillations do take place in the SN core, then the most likely detectable GWs signal should be produced over 
the time interval for which the conditions for flavor conversions to occur are kept, i. e., ~ (10 _1 — 10) ms. This 
time scale implies GWs frequencies in the band: /gws ~ [10 - 0.1] kHz, centered at 1 kHz, because of the maximum 
v production around 1 — 3 ms after core bounce (see Walker & Schramm 1987). This frequency range includes the 
optimal bandwidth for detection by ground-based observatories. For a 1 ms conversion time span the v luminosity 
reads 



AE„ 



3 x 10 51 erg = g ^ 1q5 4^E 



ATqsc 1 x 10" 3 s s 

Hence, the GWs luminosity, ^q\Ys' as a Iunc tio n of the v luminosity can be obtained from the equation 



16ttG 



dh 
~dl 



1 



= 4^ L GWs > 



(43) 



(44) 



which relates the GWs flux to the GWs amplitude, h v . In the case of GWs emission from escaping z^s this amplitude 
is computed from the expression (Epstein 1978; Burrows & Hayes 1996; Miiller & Janka 1997) 



hij 



2G_ 

c 4 R 



L 



t-R/c 



dt' L v {t') a(t') e l ®e 



3 ' 



(45) 



where ® ej represents the GWs polarization tensor. 

To attain order of magnitude estimates one can transform Eq.(45) to get the amplitude of the GWs burst produced 
by the non-spherical outgoing front of oscillation-produced v s as (Mosquera Cuesta 2000; Burrows & Hayes 1996; 
Miiller & Janka 1997) 



2G 



K 
K ~ \A\ 



A (At x L u xa) 



(46) 



55kpc 



R 



AT 



3 x 10 54 ^£ 



\0.lJ 



where \A\ = 4 x 10 23 Hz 1//2 is the amplitude. 4 Equivalently, that amplitude can be reparametrized as 



K ~ \A\ 



55kpc 



R 



AE, 



GWs 



10- 7 M Q c 2 



1/2 



10 d Hz 

/GWs 



1/2 



(47) 



A GWs signal this strong will likely be detected by the first generation of GWs interferometers as LIGO, VIRGO, 
etc. Its imprint in the GWs waveform may resemble a spike of high amplitude and timewidth of <~ ms followed by a 
Christodoulou's memory (Mosquera Cuesta 2002). From Eq.(47) the GWs luminosity turns out to be 



J GWs 



10 



50-48 er S 



3 x 10 54 ^£ 



a 



0.1 - 0.01 



while the GWs energy radiated in the process yields 

£ GWs = L GWs x AT osc 



10 



47-45 



erg 



(48) 



(49) 



This is about 10 5 3 x L^ iS - the luminosity from v diffusion inside the PNS, as estimated earlier (Mosquera Cuesta 
2000; 2002). 

EDITOR PLEASE PLACE TABLE 5 AROUND HERE !!! 



4 One must realize at this point that the high value of the anisotropy parameter here used is consistently supported by the 
discussion regarding the neutrino coupling to rotation and magnetic field presented above. 
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VI. CONCLUSION 



One can see that if v flavor conversions indeed take place during SN core bounce inasmuch as they take place in 
our Sun and Earth (Smirnov 2002), then GWs should be released during the transition. The GWs signal from the 
process is expected to irradiate much more energy than current mechanisms figured out to drive the NS dynamics at 
birth do. A luminosity this large (Eq.(44)) would turn these bursts the strongest GWs signal to be detected from 
any SN that may come to occur, futurely, on distances up to the VIRGO cluster, R ~ [10 — 20] Mpc. It is stressed 
that this signal will still be the stronger one from a given SN, even in the worst case in which the probability of v 
conversion is three orders of magnitude smaller then the estimated in the present paper. In proviso, we argue that 
a GWs signal that strong could have been detected during SN1987a from the Tarantula Nebula in the Milky Way's 
satellite galaxy Large Magellanic Cloud, despite of the low sensitivity of the detectors at the epoch. In such a case, 
the GWs burst must have been correlated in time with the earliest arriving neutrino burst constituted of some active 
species given off during the very early oscillation transient where some v e s went into v^s, v t s or ^ s s. Thenceforth, it 
could be of worth to reanalyze the data collected for from that event taking careful follow up of their arrival times, if 
appropriate timing was available at that moment. 
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EDITOR: THIS IS FIGURE 1 !!!! 
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FIG. 1. The mechanism for v conversions. Firstly, the state of a system, vs in the PNS, is defined by the energy Ea for the 
time interval ts — tA- After the flavor transitions that system is represented by a reduced energy Eb over time to — tc- The 
total energy transferred to the new species is AE^ a ^ s = Ea — Eb, and thus the v luminosity over the transient is given by 

T _ E A -E B 
Uv ~ *C-*B ■ 



1G 
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EDITOR: THIS IS TABLE 1 !!!! 

TABLE I. Time scales (in ms) called for in this paper. 



Supernova v Neutrino Supernova Supernova 

Thermalizat. Oscillations Deleptonizat. Core-Bounce 

A rrithcrmal A rpU A rriPNS A mPNS 

<-±J-v ai osc '-i-'delept ^'bounce 

~ 500 ~ 3 - 10 ~ 10 ~ 20 
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TABLE II. Parameters of the PNS used for in this paper. 



Radius 

Rpns 
~ 80 km 



Mass 

MpNS 

' 1.3 M 



Density 

Ppns 
3 x 10 n g cm" 3 
2 x 10 14 g cm~ 3 



//-sphere 

^PNS 

35 - 60 km 



Ang. Mom 
a 

~ 0.9 



1 



EDITOR: THIS IS TABLE 2 !!!! 
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ccccc 

TABLE III. Matter densities relevant for two-u oscillations. Notice that — > v T oscillations take place as in vacuum. The 
coupling Fermi constant is: G% = 5.29 x 1CP 44 cm 2 MeV~ 2 , or equivalently Gf = 1CP 49 erg cm -3 . 







V e — > V a 


— > v T 




A 

2V2G F E 


N e 


N e - \N n 





-iiV 

2 n 



EDITOR: THIS IS TABLE 3 !!!! 
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TABLE IV. Parameter range for GWs emission from vs oscillations in SNe. The symbol DM stands for v a as a Dark Matter 
candidate. 

I'-sphere v Luminosity Oscil. Length v Velocity Am 2 (u a <-> v a ) Am 2 (c a <-» v s ) Am 2 (v a <-» v s ) 

Rpns L(v-Diff.) y oBC V v -Diff. (SNO,KamLAND) (WMAP) (v a DM) 

~ 35 km ~ 10 53 erg s" 1 ~lkm ~ 10 9 cm s" 1 ~ 10 -4 eV 2 - 1 eV 2 ~ 10 4 eV 2 



EDITOR: THIS IS TABLE 4 !!!! 
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TABLE V. GWs estimates from v oscillations and diffusion. 



GWs Freq. 


GWs Energy 


L v Asymm. 


GWs Energy 


/gw b 


£gw(^:) 


a 


E GW (v - Diff.) 


[10 - 0.1] kHz 


[ 1Q -6 _ 1Q -8] Mq c 2 


[0.1 -0.01] 


[1Q- 10 - 1Q- 13 } M c 2 



EDITOR: THIS IS TABLE 5 !!!! 
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